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The vertex-cover problem on the Hanoi networks HNS and HNS is analyzed with an exact renor- 
mahzation group and parallel-tempering Monte Carlo simulations. The grand canonical partition 
function of the equivalent hard-core repulsive lattice-gas problem is recast first as an Ising-like canon- 
ical partition function, which allows for a closed set of renormalization-group equations. The flow 
of these equations is analyzed for the limit of infinite chemical potential, at which the vertex-cover 
problem is attained. The relevant fixed point and its neighborhood are analyzed, and non-trivial 
results are obtained both, for the coverage as well as for the ground-state entropy density, which 
indicates the complex structure of the solution space. Using special hierarchy-dependent operators 
in the renormalization group and Monte Carlo simulations, structural details of optimal configura- 
tions are revealed. These studies indicate that the optimal coverages (or packings) are not related 
by a simple symmetry. Using a clustering analysis of the solutions obtained in the Monte Carlo 
simulations, a complex solution space structure is revealed for each system size. Nevertheless, in the 

■ thermodynamic limit, the solution landscape is dominated by one huge set of very similar solutions. 

t3 ■ 

^ i. introduction 

We study the vertex-cover problem [U, Q on the recently introduced set of Hanoi networks (3l-[5|(42|. An optimal 
' , ■ vertex cover attempts to find the smallest set of vertices in a graph such that every edge in the graph connects to 
^ at least one vertex in that set. It is one of the classical NP-hard combinatorial optimization problems discussed in 
Q , Ref. [6J . The problem is equivalent to a hard-core lattice gas [7] , in which any pair of particles must be separated by 
' at least an empty lattice site. The vertex-cover problem has recently attracted much attention in physics, because in 
ensembles of Erdos-Reny random networks [sj, phase transitions in the structure of the solution landscape were found 

■ that coincide with a polynomial-exponential change of the running time of exact algorithms [1, 2]. 

^ I During the past decade, alternative ensembles of random networks have attracted the attention of physicists. Well- 
i known examples are Watts-Strogatz small- world networks [9(] and scale- free networks (lol - [l3i . These networks exhibit 
' more structure and describe the behavior of real networks much better than Erdos-Reny networks Also, physical 
C^^ , systems (such as the Ising model fisl , [l6j ) that exist on these more complex network or lattice structures behave 

' differently compared to regular (hyper-cubic) lattices or random networks. 
CO , Hanoi networks mimic the behavior of small-world systems without the usual disorder inherent in the construction 
' of such networks. Instead, they attain these properties in a recursive, hierarchical manner that lends itself to exact 
real-space renormalization [17]. These networks do not possess a scale-free degree distribution; they are, like the 
' , [ original small worlds, of regular degree or have an exponential degree distribution. These Hanoi networks have a more 
J> physically desirable geometry [18], with a mix of small- world links and a nearest- neighbor backbone characteristic of 
lattice-based models j^. 

^ . For the vertex-cover problem considered here, or the equivalent hard-core lattice gas, it is difficult to find metric 

, structures with a non-trivial solution. For instance, hyper-cubic lattices are bipartite graphs that always have an 
obvious unique and trivial solution without any conflicts. Of the planar lattices, the triangular one is certain to 
exhibit imperfect solutions (i.e., there will be edges requiring multiple coverings for any solution), but any such 
solution is translationally invariant and can be easily enumerated, leading to a vanishing entropy density. Similarly, 
a fractal lattice such as the Sierpinski gasket, say, only has trivial solutions of that sort. Both of these examples are 
given in Fig. [TJ In contrast, we find an extensive ground-state entropy here, similar to the anti-ferromagnet on a 
triangular lattice [19]. Yet, our ground states do not appear to be the result of any symmetry relation. Thus, the 
study of the vertex-cover problem on the Hanoi networks affords simple, analytically tractable examples of coverages 
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Figure 1: (Color online) Vertex covering for (a) a triangular lattice and (b) a Sierpinski gasket. In both cases, the optimal 
coverage (large dots) is imperfect (i.e., some edges possess double coverings). Yet, these solutions either are unique, as for the 
Sierpinski gasket, or possess a finite symmetry, such as the possible translations on the triangular lattice, both cases leading to 
a vanishing entropy density. For both lattices it is easily seen that the asymptotic coverage is |. In the case of the triangular 
lattice, the unit cell (dashed red box) contains two vertices completely and shares half of eight vertices with other cells, i.e., it 
has effectively 2 + § = 6 vertices of which 1 + | = 4 are covered. The unit cell in the Sierpinski gasket contains 3 + | vertices 
of which the three fully contained ones must be covered. 

that have nontrivial entropy densities. In fact, analytically we found merely an approximate algorithm to generate 
(and enumerate) the set of all solutions whose true cardinality we can determine at any finite system size only by 
exact renormalization. 

Using branch-and-bound algorithms, we enumerate exact solutions however, due to the exponentially growing 
running time of this exact algorithm, we are restricted to rather small system sizes. Hence, for most of the numerical 
studies performed here, we use Monte Carlo simulations f20| to generate the solutions and clustering algorithms to 
elucidate their correlations (2l] |. 

Previous work Q has focused on averaged properties on locally tree-like (mean-field) networks using the replica 
method, unearthing interesting phase transitions for the problem. Thus far, there are only a few investigations into 
the statistical mechanics of the vertex-cover problem on more complex networks. In a study of randomly connected 
tetrahedra [22], glassy behavior was observed. When introducing degree-correlations, it was found that the vertex- 
cover problems becomes numerically harder ^23.] . 

This paper is organized as follows. In Sec. |ll]we review the properties of the Hanoi networks. In Sec. IIIII we briefly 
recount the relevant theory for a thermodynamic study of vertex cover in terms of a hard-core lattice gas. In Sec. IIVI 
we develop the renormalization group treatment of the lattice gas, with most of the technical details deferred to the 
Appeirdix IVIl and its applicatioir to the Hairoi iretworks HNS and HNS. A detailed numerical study of the problem 
follows in Sec. |Vl We present our coirclusions and an outlook for future work in Sec. IVIl 

II. GEOMETRY OF THE HANOI NETWORKS 

Each of the Hanoi networks possesses a simple geometric backbone, a one-dimensional line of sites < n < N — 
2*^-1-1 [1, Most importantly, all sites are connected to their nearest neighbors, eirsuring the existeirce of the 
Id-backbone. To generate the small-world hierarchy in these networks, consider parameterizing airy integer n (except 
for zero) uniquely in terms of two other integers (i, j), « > 1, via 

n = T~'{2j + 1), (1) 

where i denotes the level in the hierarchy and j labels consecutive sites within each hierarchy. For instance, i = 1 
refers to all odd integers, i — 2 to all integers once divisible by 2 (i.e., 2, 6, 10,...), and so on. In these networks, 
aside from the backbone, each site is also connected with some of its neighbors within the hierarchy. For example, we 
obtain a 3-regular network HN3 (best done on a semi-infinite line) by connecting first the backbone, then 1 to 3, 5 



Figure 2: Depiction of the 3-regular network HNS on a semi-infinite line. Note that HN3 is planar. 



to 7, 9 to 11, etc, for i = 1, next 2 to 6, 10 to 14, etc, for i = 2, and 4 to 12, 20 to 28, etc, for i = 3, and so on, as 
depicted in Fig. [51 Previously (sj , it was found that the average chemical path between sites on HNS scales as 



iHN3 



Vl (2) 



with the distance / along the backbone. 

While HN3 is of a fixed, finite degree, there exist generalizations of HNS that lead to new, revealing insights into 
small- world phenomena 0, 0, H^l • For instance, we can extend HNS in the following manner to obtain a network of 
average degree 5, hence called HN5. In addition to the edges in HNS, in HNS we also connect each site in level i {i > 2, 
i.e., all even sites), to (higher- level) sites that are 2*~^ sites away in both directions. Note that Eq. (jlj implies that 
the nearest neighbors of a site i within its hierarchy are separated by a distance of 2 x 2*~^. The resulting network 
HNS remains planar but now sites have a hierarchy-dependent degree, as shown in Fig. [3] To obtain the average 
degree, we observe that 1/2 of all sites have degree S, 1/4 have degree S, 1/8 have degree 7, and so on, leading to an 
exponentially falling degree distribution oi V {a = 2i + 1} (x 2^'. Then, the total number of edges L in a system of 
size iV = 2*^ -f 1 as shown in Fig. [3] is 



fe-i 



2L = 2 (2/c + 1) + {2i + 1) 2''-' = 5x2'' -i, (3) 

i=l 

where the expression outside the sum refers to the special case of those three vertices at the highest levels, k — 1 and 
k. Any other choice of boundary conditions may vary the offset in Eq. ([3]), but not the average degree, which is 

o r 

(a) = — ^ 5. (4) 

In HNS, the end-to-end distance is trivially 1 (see Fig. [3]). Therefore, we define as the diameter the largest of the 
shortest paths possible between any two sites, which are typically odd-index sites farthest away from long-distance 
edges. For the = 33 site network depicted in Fig. [31 for instance, that diameter is S, measured between sites 3 and 
19 (starting with n = as the left-most site), although there are many other such pairs. It is easy to show recursively 
that this diameter grows as 

d"N5 ^ 2 [k/2\ + 1 - fog2 N. (S) 



Other variants of the Hanoi networks are conceivable. For instance, a non-planar version has been designed [25|, [26 
but that network happens to possess only a unique, alternating covering of i and is not considered here. 



III. VERTEX-COVER PROBLEM AS A HARD-CORE LATTICE GAS 

Vertex cover is a well-known NP-hard combinatorial problem [23, that consists of finding a minimal covering 
of the vertices of a network in such a way that each edge is covered at least once. Formally, for a graph G = (V, E), 
with V being the set of vertices and E C F'-^-' the set of edges, a vertex cover V' is a subset of V with the property 
that for each (undirected) edge {i, j} G E cither i G 1/' or j G V' . A minimum vertex cover Vinin is a vertex cover of 
minimum cardinality |Vmin|- 



Figure 3: (Color online) Depiction of the planar network HNS, consisting of an HNS core (black lines) with the addition of 
farther-reaching long-range edges (shaded lines). Note that sites on the lowest level of the hierarchy have degree 3, then degree 
5, 7, etc, comprising a fraction of 1/2, 1/4, 1/8, etc., of all sites, which makes for an average degree of 5 in this network. (There 
is no distinction made between black and shaded lines in our studies here.) 



As shown in Ref. [7], the vertex-cover problem can be formulated alternatively as a hard-core repulsive lattice gas 
problem. In this formulation, the uncovered vertices of the covering problems correspond to the actual gas particles. 
These particles have a hard-core repulsion such that they can not occupy neighboring lattice sites, i.e., they cannot 
simultaneously vie for the same edge. Interpreting these particles as the voids of the covering problem implies that 
no edge may be left uncovered on both ends. Accordingly, all properties of the minimum cover problem derive from 
the ground state of the lattice gas at its highest packing. 

The grand canonical partition function for such a lattice gas is generically given by 

= ■■■ exp I ^^a;J (1- XiXj), (6) 

a:o = {0,l} 2;iv = {0,l} \ i=l / (i.j) 

where the product extends over all edges of the graph and exerts the hard-core repulsive constraint. The chemical 
potential ^ is provided to regulate the density as gas particles get packed into the system. Since maximal density of 
the gas implies minimal coverage of all edges, we are looking for the configurations in the limit /i — > oo of the gas. 
The quantities we seek are the thermodynamic limit {N — >■ oo) of the packing fraction for the lattice gas, 

"M - ;^(|:-.) 4|;1"HM, (7) 



and the entropy density of such configurations, 

siHf^)) = ^(l-f^-^ynEi^^). (8) 

It has also been shown in Ref. [t] that one can extract the corresponding properties of the minimal vertex coverage 
from these in the /i — ^ cx) limit. For the coverage density, this corresponds simply to the void density of the gas, 

Cmin = 1- lim V (fl) , (9) 

>oo 

and the entropy density of optimal coverages is simply equal to that for the lattice gas: 

SVC (Cmin) = S {u ^ 1 - Cmin) ■ (10) 

Due to the hierarchical structure of the Hanoi networks, we will also introduce level-specific chemical potentials /i^, 
for example, to extract information about the coverage with respect to the level of the hierarchy (i.e., the range its 
small-world edge attains) that a vertex may reside in. The corresponding derivations are presented in the Appendix. 
Throughout, we will find it often convenient to express the chemical potentials as an activity variable. 



(l<z<fc), (11) 



such that /ii — )■ oo corresponds to the somewhat more tractable limit — ^ 0. 
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IV. RG FOR THE HARD-CORE LATTICE GAS ON HANOI NETWORKS 



The renormalization group (RG) as applied to the lattice-gas problem developed here contains a few unfamiliar 
features. Thus, we have to elaborate to a significant extend on the procedure. Although ultimately the RG will heavily 
rely on procedures used for Ising spin models, initially we will have to rewrite the grand canonical partition function 
of the lattice gas in an appropriate form. To this end, the purpose of the first step of the RG (already eliminating 
half of all sites) is to generate the initial conditions for the subsequent canonical partition function analysis, in which 
the usual coupling variables depend in a complicated way on the chemical potential ji instead of a temperature, and 
the apparent "spin" variables are in fact Boolean, Xi € {0, 1}. 

We have to rewrite the generic partition function in Eq. ^ for the special case of the Hanoi networks. To access 
more details of the solutions, we will take the opportunity to generalize to the case of a hierarchy-specific chemical 
potential /ii for 1 < i < k, where = 2'^ -I- 1 is the size of the system. (For the RG, it is natural to consider the Hanoi 
network with an open boundary both at node and at node 2^] for a system with periodic boundaries on a loop, both 
of these nodes would become identical and N = 2*' would be the size of the system. Of course, either choice results 
in identical thermodynamic averages.) 

First, we rewrite the hard-core repulsive factor in Eq. ^ as separate products, one for the long-range edges and 
the other for the backbone edges. 



K 2'" + i"' 

n n (1 

i—l n—1 



'^2'-i(n-l)2;2»- 



fe-1 2"' 



i=l 1 = 1 



(12) 



The case K, — 1 corresponds to HNS, with a simple, one-dimensional line of edges connecting all sites in the backbone 
sequentially. In turn, for HN5 we set /C = fc, with each i > 1 referring to the layers of those edges that connect along 
the backbone only every second site, every fourth site, every eight site, etc., as shown in Fig. [31 Note that in Eq. ([Tz 
we have used the decomposition of the sites in the network implied by the renumbering in Eq. ([1} . 
By the same token, we re-order the summation in Eq. (|6]) as 



t{0) 



i{N) ' 



(13) 



Xq 



E 



-Xg-X^^l -X k 



xo,x^k-i ,X^k 



fe-1 2"- 



n n E E 



'-2'-i(4i-3) ■^2'-i(4!-l) 



i=l 1 = 1 X„i-i 



2*-^(4i-3) ■^2^--^ (41-1) 



where we have simplified the notation on the sums to mean = X]xe{o i}' course, Eq. (jl3p has to be understood 
in an operator sense, i.e., the summations extend to all site- variables that match the indicated index. Here, we have 
also allowed for a site-specific chemical potential. It is our goal to extract local packing information, not for each 
site, but for all vertices within a specific hierarchy, where i{n) refers to the chemical potential in the z-th level that 
the vertex n is associated with according to Eq. ([T]). Naturally, the sites at the highest level k of the hierarchy 
{xo,X2k-i,X2k) require special consideration. 

In this parametrization of the indices, the products in Eq. can be combined with those of the second factor in 
Eq. (fT2l) . Both refer to the small- world edges in all levels of the hierarchy and are naturally expressed in a hierarchy- 
conform manner. Hence, we find for the grand-canonical partition function defined in Eq. ([S]) on a Hanoi network 
with k levels in the hierarchy: 



(mi 



,mk) 



E 



TO, 



— XQ—x^k-l —x^k 



5/c (m2, . . . ,mfc_i) Y[ Q {mi,X2(2j-2),X2{2j-i),X2{2j)) , (14) 



where we have defined the operator for the weighted summation on HN3 and HNS, respectively. 



5] 



HN3 



n n E E 



TO, 



-x^,-i 



2'-i(4i-3) ^2'-l(4i-l) /-I \ 

(1 - a^2'-i(4i-3)a^2'-i(4i-l)j , 

2'-i(4!-3)^'^2'-i(4i-l) /-■ \ 

(1 - a;2i-i(4i-3)a;2i-i(4i-i)j 



fc-l 2" 



5, 



HN5 



n n E E 



TO,- 



-x^,-i 



^2*-^(4;-3) ^2^-l(4Z-l) 



1=2 1 = 1 x^i-i^ 

X (l - a;2'-i(4(-4)a^2'-i(4i-3)) (l - a^2'-i(4i-3)a^2'-i(4(-2)) (l - a^2'-i(4(-2)a^2'-i(4(-l)) (l - a^2'-i(4i-i)a^2'-i(4i)) 
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Note that these operators only sum over ah even-indexed variables (i.e., i > 2). To obtain a renormalizable form for 
the partition function it is necessary to trace over the lowest level i = 1 of the hierarchy, i.e., to eliminate all odd-index 
variables. For both, HNS and HN5, this results in an identical structure, defined as 

(Mi,a;2(2j-2),a;2(2j-i),a;2(2j)) = ^ ^ m'''-''-^''^"-^ {1 - X4j-3XAj-i) (15) 

X4j — S X4j — i 

(1 - X^j^iX^j^s) (1 - X4j^3X4j^2) (1 - X4j^2X4j-l) (1 " XAj-iX4j) , 

= 1 + (l - a;2(2j-i)) (2 - 0:2(2^-2) - 2:2(2^)) • 

In Appendix A, we show how to recast O in an Ising-like form with a sufficient number of renormalizable parameters. 
We can simplify the grand partition function in Eq. (|14p further by combining the products and writing 



2^=) (mi,...,TOfc) = J2 



TO I. 



-Xo — X^k-l —X2k 



fe-2 2" 



n n E E 



1=2 1 = 1 X 



2^(42-3) ■^2^(4i-l)_ 



n ^1 {^4(21-2), X4{2l-l),X4(2l)){'^G) 



1=1 



where the explicit expression for C| is also derived in Appendix A for both, HNS and HNS, which allows us to drop the 
subscript label. In either case, the RG recursion equations now result from imposing the recursive relation between 
hierarchies, 

Ci+i (a^2'+i(2i-2),2;2i+i(2i-i),a;2i+i(2;)) (17) 

= E E C?^^ {^2' (4l-A),X2i (41-3)^X2' {41-2)) {X2i (41-2), (41-1), (41)) 1 

^2^(41-3) ^2^(4Z-1) 

which are derived in Appendix A. There, Figs. [T5] and 1161 also provide a graphical representation of Eq. (|17p . 



A. Analysis of the RG Recursions 

We find that the RG recursions that follow from the previous discussion, which are given explicitly in Eqs. ((52|) for 
HNS and in Eqs. ((54|) for HN5 for the hard-core lattice gas model, have only two trivial fixed points. There is a stable 
low-density fixed point for all /x < 00, i.e., to > 0, and an unstable fixed point at full-packing for fi = 00, i.e., to = 0. 
Note that in this part of the analysis we are concerned with global properties, and thus, ignore differences between 
the hierarchical level by setting = to throughout. 



1. Analysis for HN3 

The limit m — > of the recursions in Eqs. ([52l) for initial conditions given in Eqs. ([50|) is difficult to handle. Except 
for Ki, all other parameters are either diverging or vanishing in Eqs. (|18ll for that limit. To achieve a clearer picture, 
we evolve the recursions once and obtain 

24 8 ^ to2 15 , 25 ^ 4 

m — , 12^ 17, C2 ^ 1^2- -5-, A2 - — , A2 - -— . 18 

5 3 8 8 24 25to 

In fact, further revolutions in the recursions seems to preserve this picture: Ci scales with a rapidly growing power 
of TO, while all other parameters and — toA^ become finite for to = at any order i. Thus, we replace A with A 
and subsequently set to —> in Eqs. ((52|) yielding 

'^i+i ^ — — ' li+i ^ liVii^i, Vi+1 ^ — -2' 

^ (1 + K,) 

. (1 + ^^.) , (1 + ^.)' X 2^fA, 

(2 + 7iKfA,) (1 + K,) 

At its core, the two recursions for k, and A have become independent of all the others. The m — Q fixed-point itself is 
then dominated solely by the stationary solution of their recursions in Eqs. ([T^ . 

K* = -J—, \* - ^- (20) 



23-1 2t 



(21 -1) 
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Figure 4: (Color online) Plot of the value of Hi after the i-th RG-step for m = 10 ^,10 10 and 10 * (left to right). At a 
length scale ^ (m) — 2' with i = — | logj m, the behavior of Ki crosses over from the value at the unstable m = fixed point, 

K* = 1/ ^2^^^^ — 1^ ~ 1.70, to the stable m = 1 (/x = 0) fixed point at which k* = 1. 



Therefore, one finds a constant solution for ij* — 4k*/ (1 + k*)^ = 1/A* and the recursion 7i-|_i ^ ji {k* /X*) with the 
solution 7j ^ 7o2"3^ which diverges for large i. The situation for is more subtle. Numerics clearly indicates its 
decay, but this could occur consistently in two ways. First, if it were to decay such that jiAi still increases, then 
Eq. pop suggests A^+i cx l/7i, but that would render jiAi constant, which is a contradiction. Alternatively, if both 
Ai and jiAi decay, then Ai-^i ^ Ai [k* / (1 + t*)]^, yielding ~ in a consistent manner. Numerical studies 
verify that the latter solution is indeed realized. 

From the terms dropped in the to limit, we can extract a cross-over scale as follows: Achieving the limit m — > 
implies that the widely occurring term TO7i in Eqs. (|52p is considered small enough to be discarded with respect 
to terms of order unity. Hence, by identifying ^ = •\/2*('") as the correlation length within the small-world metric 
supplied by Eq. using 7i(m) ^ 1/to yields 2**^™' ~ to~2 or 

^ - eif" (21) 

as the diverging length below which the systems orders for an correspondingly diverging chemical potential, — cx). 
Indeed, for m = lO^**, for example, we find numerically that the solution veers off the unstable fixed point just below 
the i = 10th iteration; Fig. [4] demonstrates the correctness of Eq. ([21]) for any small to. 



2. Analysis for HNS 

The analysis for HNS is surprisingly subtle. Although the preceding fixed-point analysis for HNS required the 
singular limit to — as part of the consideration, after the appropriate rescaling of the parameters with to, the 
subsequent approach proceeds in a familiar fashion. HN5 obscures this approach with an additional layer of complexity, 
resulting from strong alternating effects order-by-order in the RG, as the numerics reveals. Of course, the initial 
conditions here are identical to those for HNS in Eqs. ([HI)) , with the same pathologies in the to limit. However, 
whereas those problems were essentially resolved for HNS after one RG-step and rescaling, see Eqs. (|18l) . here we find 



(22) 
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and 

16 9 1 

T^' 73 ?73 ~ 16to, Ksr^—, As--——, A3 - 16m, (23) 
16 ym 16 16to 

etc. This alternation between regular and singular behaviors of each of the parameters persists thereafter. Leaving the 
recursion for Ci aside for now, we notice that for even indices, j2n, '72n, 'mK,2m ^2m and remain finite for m — >■ 0, 
but for odd indices, this is true for m72n_i, rj^n-i/m, K2n-i, m\2n-i, and A2„-i/to. Defining 72n-i = ^7271-1, 
^2n-i = ??2ti-i/'Ti, K2n = TOK2n, ^2n-i = wA2n-i, and A2„_i = A2„-i/m, it is useful to rewrite the recursions in 
Eqs. ([54|l separately for even and odd indices. In fact, the limit m — > on its explicit appearance can now be taken 
to get 

72« = ^2n-l (2 + 72n-l), 72«-l (2 + m72(„_i)) --2772(„-l), 

2 + 72„_i _ 2 + TO72(„_i) 

m-n = l2n-\— : ^, ??2«-l = 72(n-l)- ^ 272(„_i), 

(l + 72„-ij (l + m72(„_i)) 

- T (l+72»-i)' ^ , l + m72(„_i) 1 

^2k — A2K-1 „ , _ , '«2„-i — A2(„_i)-— 7:'\2(n-l)i l^4J 

2 + 72n-i 2 + m72(„_i) 2 

. (l+72n-l)^ T 1 + m72(„_i) 1 
^^2" — ~ /f, , 7: A2n-1 — 77— r ~ , 

72«-i(2 + 72n-i) 72(«-i) (2 + m72(„_i)) 272(„_i) 

. _ 2 + 72n-i - 2 + m72(„_i) 

L^2n — l2n-l- Z 3' ^2™-! = 72(n-l)- ""J ^ ^72(n-l)- 

(l+72ri-l) (1 + »77.72(„-l)) 

Note that for the limit m — >■ we only assumed that TO72(n-i) ^ 1 for n — > oo on the right-hand set of these relations, 
which provides a correlation length from the cross-over n^o — n{m) at 72nco ~ 1/m. Eliminating all odd-index 
quantities from the equations yields 

72n = 472(„_i) (1 + ?72(n-l)j , '72n = 4772(„-l)- "3 , (25) 

(l + 2772(„-l)) 

1 + 2772(n-l) (1 + 27?2(„-l))^ _ !+%(«-!) 

'«2n — 7— T, A2„ — J— r-, l\2n — 4r72(n-l)- 



472(«-l) (1 + ??2(n-l)) 47;2(„-l) (1 + %(«-!)) (l + 2r72(„_i) 



2 ■ 



These interlacing recursions now have a simple fixed point, which derives from the only non-trivial solution of the 
self-contained ry-equation: 



73 
2 ■ 



(26) 



This implies the equally stationary value 

'\* — 

A* (l + 2?7*) 



^ ^ ^ 4r;*(l-^77*) ^ 3jV3^ ^^7) 



but we also find the asymptotically scaling 

72n ~ 70 

This provides the correlation length estimate 



(2 + 73)]" cx (28) 

V / J K2n 



i = 2"- ~ exp <^ —f- ^ \ . (29) 

^1 log2[2(2 + y3)] J 



B. Packing Fraction and Entropy 



To understand the most pertinent features of the problem, such as the optimal packing (or coverage) and its entropy, 
we have to consider the asymptotic behavior of the renormalization group parameter Ci. related to the growth of 
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the overall energy-scale, in Eq. (|19p for the initial condition in Eq. (fT5)) . Clearly, the partition function at any finite 
system size is a polynomial in e^, i.e., in powers of m^^. Both of these quantities, packing fraction and entropy, derive 
from the most divergent power in m to be found in S. To wit, we can write for m — )■ with iV = 2''' + 1, 

S(*=) - {am-")^ [l + am + bm^ + . . .] . (30) 

Then, it is S^lnS = — mSmlnS ~ Na, and we find from Eqs. (I7]and[51), 

v ~ a, 
s = In cr, 

for — > oo at TO = 0. 

Equation provides the grand canonical partition function E^''^ for 2*^ + 1 site-occupation variables in terms of 
an Ising-like canonical partition function Z'^''~^'> for only 2*^"^ -I- 1 (Boolean) spin variables. While S'^'^^ depends only 
on the hierarchical chemical potentials to^, ostensibly Z^'^^^-' depends on a tuple Ai of renormalizable couplings, see 
Eq. (|55p . in addition to any explicit dependence on to^. Of course, the couplings themselves are merely a function 
of the chemical potentials, Ai = Ai (toi), through the RG initial conditions in Eq. ([50)1 . Step by step in the RG, 
the couplings transform according to Eq. (j56p each time the system size halves, whereas the partition function stays 
invariant. Hence, we can expand on Eq. (jl6p and write 



(mi,...,mfe) = Zt'^ 1) (^11(7711), TO2,...,mfe) , 

= z'-'''^^ (^A2 (mi, TO2) , m3, . . . , TOfc^ , 



: (31) 
= Z*^) (mi,...,TOfc„i),mfe^ , 

where Z'^^-' is simply a rudimentary Hanoi network consisting of just three vertices. 

1. Results for HNS 

Specializing this discussion for HN3, we find for the rudimentary partition function Z*^^^ in this case 

(1) _ „_i -^{xo+x^k-i+x^k) -^[{xo+x^k-i'j + ix^k-i+x^k)] 

^ - ^k-1 2^ 2^^k 7fc_i 

-^{xO + X^k) -{xoX^k-l+X^k-lX^k) ^-XQX^k A-XoX^k-lX^k /qr,\ 

Vk-i ■^k-i ■ l-J^J 

For a uniform chemical potential, = to for all i, one finds that for m — >■ the partition function is dominated 
overwhelmingly by the renormalized value of C^, i.e. 

lnS('=)(^) = lnZ(i)(i'fc_i (to),to) - -lnCfc_i. (33) 
Rewriting the recursion for Ci in Eq. in this form yields 

InC.+i = 21na+ln(^) ^21na + |ln2 + ln(^), (34) 
which is easily summed up to give 

lnCfc_i = 2'=-3[inC72+ln(2TO7o)]. (35) 
With C2 ~ TO^, as listed in Eq. p^ . we get 



1 , ^(k) 1 , ^ 3. . . 1 



2fe— - 2k"— ^ - -gln(m)--ln(47o), (36) 

and comparison with Eq. (I30p produces an exact prediction for the maximal packing fraction of the lattice gas, 

3 

i/(yu ^ 00) = -, (37) 
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Figure 5: (Color online) Plot of (a) the packing fraction uvc and (b) its entropy density svc for the lattice-gas problem on 
HNS for the first few system sizes A'^ = 2*^ + 1 with fc = 2, . . . , 5 (top to bottom at m = 1) as a function of m. 




Figure 6: Plot of the relative packing per level 2Vi on HN3 for various system sizes A*' = 2* + 1 with fc = 7, 12, 17, 22, and 26, 
plotted also on a relative level-scale i/fc at m — >■ 0. Asymptotically, in large systems, all vertices in higher levels i appear to 
be just 50% packed (or covered), which is minimally necessary to cover the one small- world edge connecting such vertices. (Of 
course, each level contains half as many vertices as any preceding level and thus contributes ever less to the overall coverage.) 
This packing may well be random as such vertices are far separated between the higher levels. A significantly lower packing 
(higher coverage) is attained only at an ever small fraction of the lowest levels to account for the overall packing fraction of | 
(coverage |). 



i.e., for the minimal fraction of vertices needing cover in HNS, it is 

5 

Cmin = -r- (38) 
o 

Note that the m-dependence of C2 and of the recursion for Q in Eqs. are crucial for this result, whereas 7^ is 
independent of m and, hence, becomes irrelevant here. Unfortunately, the entropy density in turn depends not only 
on the asymptotic form for 7^ but on the iroir-trivial integration constant 70, which can not be determined from the 
asymptotic behavior of the RG flow; it is a global property of that flow and could depend on all its details. However, 
the result suggest that, at least for HNS, unlike for those lattices in Fig. [l] the entropy density does not vanish but 
attains a non-trivial value. In fact, using the recursions in Eqs. (j52p for arbitrary m and taking the m — > limit 
only in the end, we can exactly determine the constant a defined in Eq. (|30)) for the first few values of k (see Tab. HI. 
Finite-size extrapolation from the numerical evolution of the RG flow up to fc = 25 levels (i.e., system size N = 2^^) 
for a finite but small value of rn — lO^"*" predicts that svc(cmin) = 0.160426(1). (Any variation of m over 10 decades 
does not affect the extrapolation at this accuracy.) For smaller system sizes we plot the packing fraction and the 
entropy density for the entire range of the chemical potential in Fig. [5] In Appendix, we describe how to evaluate 
derivatives of the partition function, such as those leading to u and s, within the RG-scheme. There we also develop 
a method to probe the packing fraction for each level of the hierarchy; those results are plotted in Fig. [S] 

In the Appendix, we derive a partial set of recursions to approximate the number of solutions given in Tab. ID Our 
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Table I: Listing of the first few values of a and svc defined in Eqs. (f30)) and (fTO)) for HNS of size iV = 2'-' + 1. The sequence for 
the total number of optimal configurations, cr^, soon develops non-trivial prime factors. The entropy density for the coverage 
Svc only converges slowly to its numerical limit. 



k 




Svc ~ Ino" 


2 


1 





3 


7 


0.243239 


4 


37 


0.225682 


5 


718 


0.205515 


6 


193284 


0.190186 


7 


8651040480 


0.178757 


8 


11491993035377280000 


0.171438 


oo 




0.160426(1) 



Figure 7: Depiction of perfect coverings on HN3 for = 3. Of all seven solutions, we omitted the three obtained by reflection 
from these. Light-colored sites belong to the vertex cover, dark-colored sites mark particles with hard-core repulsion that 
prevents nearest-neighbor occupation. 



failure to obtain a closed set of such equations (and an asymptotic prediction) indicates the non-trivial origin of the 
entropy density. Here, we just plot the exact solutions for fc = 3 and 4 for illustration in Figs. [7] and [8l As the 
numerical results in Sec. |V] indicate, the optimal packing of the lattice gas at any finite size TV = 2*^ + 1 contains for 
any fc > 3 exactly 3x2'^^'^-|-l particles. 



2. Results for HNS 

For HN5, we find that the rudimentary partition function Z^^^^ is like that for HN3 in Eq. (|39|) . except for additional 
repulsive terms: 

XO K^fe-l x^k 

-\(xa+x^k) -(xax^k-\+x^k-\x^k^ -XQX^k /^-xax^k~\x^k /oq\ 
(1 — XQX2k-\) (1 — Xik-\X-2k) (1 — XQX-2k) . 

Hence, Eq. p3p again applies, putting the focus on the analysis of the recursion for Ci, which in its even and odd 
versions read 

2 + 72ri-l 2 -I- TO72(„_i) ^ ' 

With the results from Sec. IIV A2] at hand, when put together in the limit m 0, both recursions combine into 

C2„ - mC^(„_i) {A72„} . (41) 

The factor A'y2n, even though it grows exponentially with n, can be ignored because it does not depend on m. It is 
again easy to sum up the logarithm of this equation (for odd values of fc, in this case) to get 



^lnCfe_i - lnC2 + — Inm - Inm, (42) 



— rlnOi._i ~ — ln02H mm ^ — 

with C2 ^ rn? from Eqs. \22\ . As for Eq. (|36p . for the maximal packing fraction of hard-core gas particles, this implies 

v{li^(x>) = i, (43) 
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Figure 8: Depiction of perfect coverings on HNS for fc = 4. Of all 37 solutions, we omitted the 17 obtained by reflection from 
these. Light-colored sites belong to the vertex cover, dark-colored sites mark particles with hard-core repulsion that prevents 
nearest-neighbor occupation. 




m 
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Figure 9: Plot of the packing fraction ia/c (left) and its entropy density svc for the lattice-gas problem on HNS for the first 
few system sizes A'^ = 2* -1- 1 with k — 2, ... ,5 (with alternating behavior) as a function of m. Each entropy drops noticeably 
in the limit m ^ 0. 

i.e., for the minimal fraction of vertices needing cover in HN5, it is 

2 

Cmin = (44) 

In parallel to Sec. IIV A 11 we can obtain only the constant a defined in Eq. (|30p for the first few values of k (see 
Tab. in. By the same procedure as for HNS above, we predict here that svc(cmin) = 0.11983(1). For smaller system 
sizes we plot the packing fraction and the entropy density for the entire range of the chemical potential in Fig. |9l 
Figure [TO] illustrates the strong alternating behavior between successive levels, here in form of their relative packing 
fractions. 
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Figure 10: Plot of the relative packing per level 2Vi on HN5 for various system sizes A'^ = 2*^ + 1 with k = 7, 12, 17, 22, and 26, 
plotted also on a relative level-scale i/k at m — >■ 0. In an alternating fashion levels attain an interlaced higher or lower relative 
packing (lower or higher coverage), which varies very little between the levels and seems to converge to nontrivial values. Notice 
that the apparent closing of the gap at the highest levels results from the numerical evaluation of the RG recursions at very 
small but still finite chemical activity (here, m — 10~®) . 



Table II: Listing of the first few values of a and svc defined in Eqs. (f30)) and (fTO)) for HN5 of size iV = 2*^ + 1. The sequence 
for soon develops non-trivial prime factors. The entropy density for the coverage svc alternates and only converges slowly 
to its numerically determined limit. 



k 




Svc = In (J 


2 


2 


0.173287 


3 


7 


0.243239 


4 


6 


0.111985 


5 


159 


0.220479 


6 


1350 


0.112623 


7 


21268575 


0.131818 


oo 




0.11983(1) 



V. MONTE CARLO SIMULATIONS 



We performed Monte Carlo simulations of the lattice gas by using the grand canonical ensemble in Eq. ([6]). To 
achieve a fast convergence of the Markov chains, we used the Metropolis- Coupled Markov-Chain Monte Carlo (MCM- 
CMC) approach [29], also termed later Parallel Tempering [so] in the physics community. The idea of (MCMCMC) 
is to perform Monte Carlo simulations for n independent replicas studied at different values of the chemical potential 
fi = . . . , fin with /ii = < /i2 < . . . < /in • One allows that the replicas are exchanged via two- replica Metropolis 
steps, such that an overall detailed balance is achieved. Details of the Monte Carlo moves have been given in previous 
works, e.g. Ref. | i31i] . The parameters for the simulations performed for this work are shown in Tab. IIIll 



N 


n 


/imax tyiCS 


17 


5 


6 2 X 10** 


33 


5 


6 2 X 10** 


65 


8 


6 4 X lO" 


129 


10 


7 1 X 10^ 


257 


17 


8 1 X 10^ 


513 


21 


8 2 X 10^ 


1025 


33 


10 1 X 10^ 


2049 


53 


30 2 X 10^ 



Table III: Parameters of the MCMCMC simulations: A'^ is the system size, n is the number of different values of the chemical 
potential /j., /tmax is the maximum value of fx, and tmcs is the total number of Monte Carlo sweeps, where in each sweep each 
variable is on average allowed to flip once and n — 1 times a replica exchange is attempted. 
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Monte Carlo Simulation Results 



For comparison with the analytic calculations, we show the numerical results for the density of particles. In Fig. 1111 
the resulting largest density i/, measured at the highest value of the chemical potential fj,, is shown as a function of 
system size N for HNS and HN5, respectively. To extrapolate to an infinite system size, we have fitted [s^l the data 
to power laws of the form 

(45) 
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Figure 11: Highest density v of the lattice gas on Hanoi networks found in the Monte Carlo simulation as a function of system 
size A^. The main plot shows HNS, the inset shows HN5. The solid lines represent fits to powers laws according to Eq. (|45[l . 
see Tab. II VI The dashed horizontal line in the inset marks the value 1/3. 



The resulting values are displayed in Tab. IIVI Note that for HNS, we fitted only even powers fc, since odd powers 
result in highest densities of exactly v = ^. The resulting values Voo agree precisely with the analytical results | 
and i for HNS and HN5, respectively. Also the coefficients describing the finite-size corrections seem to be rational 
numbers ^ = | and c = — 1 for HNS and h = ^ and c — —1 for HNS. They can be understood in the following way, 
e.g., for HNS: The number of nodes is = 2^ + 1, i.e, exactly one more than a power of 2. The number of occupied 
nodes for the highest density is exactly | of the 2^ nodes plus one extra node, i.e., Nv{N) = | 2'^ + 1 = | [2^ + 1) + | 
which results in v{N) = | + | . In a similar way, the scaling for the HNS graphs can be explained, where N is 
not divisible by 3. 



Network 


Voo b 


c 


HNS 

HNS (fc even) 


0.3750000(2) 0.62500(2) 
0.333333(7) 0.3333(1) 


-1.00000(1) 
-1.0000(1) 



Table IV: Result of power law fits to the iy{N) data show in Fig. [TT] according to Eq. (|45ll . Note that for HNS, only the data 
for even powers k where used. 



Next, we go beyond the analytical calculations by studying the properties of the solution landscape via sampling 
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configurations of highest density. Hence, one must ensure that configurations exhibiting the same statistical weight in 
Eq. ([5]) are sampled with the same probability or frequency. For many systems exhibiting complex solution landscapes, 
this is quite an effort (s^ - fSGt - 

To achieve unbiased sampling here, we always stored a configuration of the highest density of a replica visiting the 
highest value /imax of the chemical potential, whenever that replica previously had visited the value ^ — in the 
(MCMCMC) scheme. It may be said that the replica has "performed a round trip". This means that before a replica 
is stored next time, it must again diffuse to /i = and return to the highest value of n [37| . Typical round-trip times 
range from around 20 for iV = 17 to around 20000 for N = 2049. To test whether this procedure yields unbiased 
sampling, we studied small systems of size N = 33, where, in principle, all solutions can be enumerated. For both 
systems, HN3 and HN5, we sampled 10^ configurations of highest density and counted how often each configuration 
was found. The resulting histograms appear very flat, see Fig. [T^l Hence the sampling seems to work very well, at 
least for Hanoi graphs. 
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Figure 12: Histogram of liow often each configuration of highest density is sampled during the (MC)^ simulation of a A'^ = 33 
node graph for HN3 (main plot) and HNS (inset). The total number of sampled configurations was 10^ in both cases. 



Next, we study the configuration- landscape of the hard-core lattice gas at the highest density. For this purpose 
we take, for each value N of the system size, a set oi K — 200 randomly sampled configurations of highest density. 
We applied a clustering algorithm to each set, to generate a hierarchical tree ("dendrogram") representation such that 
"similar" configurations are grouped closer to each other than less similar configurations. As a measure of similarity 
between two configurations {x,-"^} and {x[^^}, we simply use the normalized Hamming distance 

d({xi")},{xf)}) = l5]<5^(.,^^(„. (46) 

i 

We apply the clustering algorithm of Ward [2l|, which has already been applied to the analysis of phase-space 
structures |31, 36, 38] (see Ref. [38] for details). The resulting dendrograms are shown in Fig. [131 The configurations 
are located at the leafs of the dendrogram, at the top of each dendrogram. Arranging the configurations from left to 
right as they appear in a dendrogram, a certain order of configurations is given. Note that the order is not unique, 
since for any node of the tree, the two subtrees can be exchanged without changing the clustering. Nevertheless, 
exchanging two subtrees has no effect on the final results. Note that any set of vectors can be clustered and represented 
hierarchically in this way. This is possible even for a set of purely random binary- valued vectors. 
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HNS, N=33 HNS, N=2049 




Figure 13: Distance-distance matrices for sets of K — 200 randomly sampled highest-density configurations. The columns and 
rows are labeled by configurations; the order of the configurations in the rows and columns is the same and is obtained via 
a clustering approach (see the text). The clustering structure is visible by way of the trees ("dendrograms") which are shown 
below the matrices. The entries of each matrix are normalized Hamming distances between different configurations, shown in 
gray scale (black indicates distance 0, white indicates distance 1). 

Whether this hierarchical clustering represents the original landscape structure well, can be investigated in the 
following way. One draws the matrix of Hamming distances by using the order of the configurations to order the 
rows and columns of the matrix. If, e.g., one takes a set of suitably large, random binary- valued vectors, the resulting 
matrices would appear basically gray, showing that the order imposed by the clustering is artificial in this case. In 
Fig. [T3]the Hamming-distance matrices are shown for a couple of sample systems. For both cases, HNS and HNS, 
at small system sizes, a complex block-diagonal structures is visible, such that each visible block exhibits a similar 
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substructure. This gives the impression of a complex hierarchical organization of the configuration space. Nevertheless, 
when going to larger system sizes, the matrices exhibit much less contrast, which strongly indicates that for N ^ oo 
the solution landscape will be similar to a set of random vectors, i.e., without any complex organization. 

This result is supported when computing the cophenetic correlations, which measure the correlation between the 
Hamming distances d and the distances dc along the dendrogram 

/C= (47) 

where [. . .] is the average over pairs of configurations. Note that dc is the sum of the Hamming distances along a path 
in the tree connecting a pair configurations, respectively. 




Figure 14: Cophenetic correlations in Eq. (|47|) as a function of system size for HNS (main plot) and HNS (inset). The solid 
line displays the function ;C(iV) = 3.25iV"°-^^ 

The resulting cophenetic correlation /C as a function of system size is displayed in Fig. [M) For both cases, HNS 
and HN5, IC decreases strongly as function of system size, taking the difference between even and odd powers k for 
HNS into account. For HNS, the data is compatible with a power law IC{N) = 3.25N~'^-^^. Hence, in the limit of 
infinite system sizes, the hierarchical structure imposed by the clustering is not correlated to the actual Hamming 
distances. This shows that the landscape of highest-density configurations appears to be simple for both HNS and 
HN5, in strong contrast to the vertex-cover or lattice-gas problem on random graphs |Sli) . 

VI. CONCLUSIONS 

We have succeeded in obtaining the optimal vertex coverage or packing fraction for the Hanoi networks HNS and 
HNS using the renormalization group. Our Monte Carlo simulations allowed us to confirm those results and extend 
them to any finite size. We have also obtained the entropy to arbitrary accuracy. We have shown that it is extensive 
and likely non-trivial in the sense that there is no simple generator to provide the set of all optimal configurations, a 
remarkable result for such a simple, planar network. It is even more remarkable that for each given size the set of all 
possible solutions has a complex hierarchical structure, as visible from clustering the states and considering distance- 
distance matrices. Nevertheless, an analysis of the cophenetic correlations shows that in the thermodynamic limit. 
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a set of random-vector-like solutions dominates entropically and makes the solution landscape thermodynamically 
simple. 

While there are no phase transitions in this problem, the Hanoi networks would allow one to study analytically an 
interesting percolation transition when considering an interpolation between the network's one-dimensional backbone 
alone (a simple bipartite lattice with just two perfect solutions of 1/2 coverage) and the full network (with an extensive 
set of frustrated optimal solutions of coverage 5/8 for HNS or 2/3 for HN5) by adding the small- world edges with a 
probability p. As a technical achievement, we derived the renormalization group equations for hierarchy-dependent 
observables to obtain, for instance, the packing fractions provided by each level of the hierarchy in the network. 
Here, these observables merely reveal that higher levels of the hierarchy become very uniform (even if alternating) 
in coverage, while most of the interesting structure resides with the majority of variables at a few lowest levels, in 
accordance with the numerical study of the ultrametric relation between solutions. Nevertheless, similar techniques 
might be useful to pro vide insights into the "patchy" nature of ordering on whole classes of hierarchical networks in 
other problems [H, [H, [H, [H E3 ■ 
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Appendix 

A: Determining the RG-Recursion Equations 



In the derivation of the recursive form of the partition function in Sec. IIV[ we use Eq. ([T5|) to transform O into the 
Ising-like form with Boolean variables x, y, z 

e(/xi,a;,y,z) = 1 + e^^l - y) (2 - x - z) . 

= exp [21 + \p \(x + y) + (y + z)\ + iff {x + z) + K [xy ^ yz) ^ Lxz + Dxyz 

= c-27-^[("+^'+(^+^\-^("+")«;-(^^+^^)A^"^A^^^ (48) 
where we have defined the convenient "activity" parameters 

C^e-', 7 = e-^, v^e'", n = , A = e"-^, A = e"^. (49) 

Equation ((i5|) matches Eq. ([T5|l for the choice of 

mi 2 + mi 7111(2 + mi) 
2 + nT-i mi (1 + mi) 

1 + mi (1 + TOi)^ mi (2 + mi) 

2 + mi mi (2 + mi) [\ rai) 

(with mi — e~^^), which serves as the initial conditions for the renormalization-group flow for both HNS and HNS. 

In terms of these renormalization-group parameters one can then show for HNS that the "sectional" partition 
functions have to be written as 



2 -a-b -^l{x+a) + ia+y) + (y+b) + {b+z)] - ^[ix+y) + {v+z)] 

a b 

[xa+ay+yb+bz] —(xy+yz) * —{xay+ybz) 



>^-yxy^yzj^-y.ay^yuzj _ ^ ^^^^ 

^-1 -i[i'^+V) + (y+^)]~i(.'^+^)^-ixy+Vz) ^-xz /^-xyz 
'^i+lli+l 'li+1 ! 
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Figure 15: Depiction of the graph-lets associated with the sectional partition function i^l in Eq. (|51|l during one RG step on 
HNS. The step consists of tracing out odd-labeled variables Xn±i (taking into account the hard-core constraint relevant at this 
level) in (a) and expressing the renormalized couplings (7' , r]' , k' , A', A') , in (b) in terms of the old couplings (7, 77, k, A, A). To 
save space, the one-point couplings ("bond magnetizations" [41| ') 7 and r) have been omitted. These drawings summarize the 
calculations in Eqs. (|51|l and ((52]). 



for which we have depicted the tracing operation graphically in Fig. [151 This operation requires that, for HNS, the 
renormahzed quantities at i -|- 1 be expressed in terms of those at i with the RG recursions 



2 -f TOi+i7i 
(2 + mj+i^j) (2 + mi+i7,K,) 
{1 + Ki+ mi+i^^Kif 

(1 + K, + mj+i7,Ki)^ 
Ki (2 + mi+i7,:) (2 + mj+i7^Kj) 

For HN5, we obtain, correspondingly. 



Ai-i-i 



Ki+i — A^Ai 



2 + nii+ijiKi 
(2 -I- mt+i7tKt) {1 + Ki+ mj+ijiKj) 
(2 -I- mi+i7j) (1 -I- KiAi + m^+ijiKfA^) ' 



(52) 



(2 + TOt+i7t) (1 -)- K.Aj + TOj+i7,Kf A,)^ 
A; (2 + mi+i-fiKf Ai) {1 + Ki + ■m.,+ijiKiY 



ri(^ „ ^\ - r-'^„,-a-b-h\-i^+a)+{'^+y)+{y+b)+{b+^)]-h[(x+v)+{y+^)\ 



(53) 



a b 



„ ri-i -i[(^+v)+iy+^)]-i(^+^) ^-(xy+y^) \-xz >.-xyz 
~ "-"i+i'i+i 'li+i %+i '^i+i ^i+i 1 

a procedure that is graphically depicted in Fig. 1161 Those extra repulsion terms in HN5 then lead to dramatically 
simpler RG recursions than in Eq. ([52]): 



rrii+ijiCf 2 + mi+i7i 



2 -I- mi+i7j 
(1 + nii+ijiKi) 

Ki+l — Ai ^ Ai+l 



TO 



rrii+i^i (2 + TOi+i7i) 



\2 ' 



(1 -I- TO^+l7^) 



Ai+i 



(1 + TOi+i7j) 

TOi+i7i (2 -I- TOi+i7j) 



(54) 



(2-|-mj+i7j) ' TOj+i7i (2-|-TOi+i7i)' " (1 TOi+i7i)^ 

For the discussion in Appendix B, it is useful to defined the vector of renornializable parameters, 

(toi, . . . ,TOi) = (Ci,7i,77i,Kj, Ai, Ai) , (55) 

where at each level of the RG i depends implicitly, through the renormalized parameters, on the first i values of the 
chemical potentials, as in Eq. ([TS]) for the initial case i = 1, for example. In the analysis, we will symbolically refer to 
these renormalization group equations formally as a (nonlinear) operator, 



(toi, . . . ,TOi+i) = Tlrn^ + i (toi , . . . , TO^ ) 



(56) 



highlighting the fact that the RG transforms depend explicitly on the parameters to^+i. 
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Figure 16: Depiction of the (exact) RG step on HN5. This step is identical to that for HNS in Fig. [15] aside from the additional 
hard-core repulsive terms in (a) between Xn±2 and a;„ which is relevant for the current RG step, and in (b) between Xn-i and 
a;„+2 which contributes at the next level of the RG. 



B. Hierarchical Packing Fraction 



For later use, we follow convention in defining the Jacobian matrix derived from a formal derivation of the renor- 
malization group equations as defined in Eqs. (|55|56p . 



A; = 



9li 



d {C^+x , 7,+i , , , Aj+i , A^+i ) 



(57) 



Using the fundamental statement for the grand partition function S*-*^) of the unrenormalized system (or the free 
energy /f*^' = 2"*^ In S^'^' , instead) in terms of the renormalized partition functions Z'*^'^) in Eq. (|3ip , we can find for 
the specific packing fraction in the z-th level of the hierarchy 



)fc ( X! 



2k 



a;2'(2j-i) / = 



dm,, 



InSW, 



implicitly defining the hierarchy-specific chemical potential nii = e^* in the form of the vector 

rh = {mi, 7712, • ■ • , ruk) ■ 
Applying such a derivative to the sequence in Eq. ((3T|). we obtain for 1 < z < A; 



- — InS^'') (mi,m2, 
am, 



dm, 



Ak^i (mi, . . . ,mfc_i) ,mfc 



51nZ(i 



Ak-i,mk 



dA, 



k-l 



(58) 



(59) 



o— . (60) 

ami 



We can understand the progression of derivatives in Eq. (pO)) from the result in Eq. ([5 

ST? 

dAi ' 



dm, 



—n 

dm, 



(mi, . . . ,m/_i) 



,0, 



A_i (mi, . . . ,mi_i) 



Ai-i o 



dAi-i(mi,...,mi-i) 



i = I, 
i < I, 

i > I, 



(61) 



using, from Eq. (|57p. the matrix 



A, 



dn. 



mi + l 



dA 



Ai (mi, . . . ,m/) 



(62) 



Note that the distinction between the implicit and explicit derivatives in Eq. ([61]) results from the explicit occurrence 
of mi just that once in the i-th RG step in the recursions and that afterward the parameters being renormalized 
depend implicitly on m^. Thus, application of the relation in Eq. (|6ip . repeatedly for alH > i and once, finally, for 
I = i, yields 



- — In.^^'') (mi,...,mfe) 
dm, 



ainz(i) 

dA 



(i'fe_i,mfc) o 



A 



k-2 ° 



A, 



dm 



A_i (mi, . . . ,mi_i) . (63) 



21 



Now it is easy to set all chemical activities equal, TOj = m. with 1 < i < fc, irrespective of which hierarchy was 
targeted, to get 



- — In^'^''^ (mi, . . . ,mk) 
ami 





r 91nZ(i) 






dA 






Mfc-i,mJ 




dm 



Ak-i,m ]o^{ Ak-2 ]o...oW {AAo^{ ) , l<i<k 



dm 



^ k. 



(64) 

We can relate this procedure back to that for the total occupation defined in Eq. ([7]) using a uniform m. To this 
end, we define an extended vector of parameters with an explicit m-dependence 



A\ = [Ai,m]^{Ci,-/i,r].i,Ki,Xi,Ai,m). 



Then 



dm 



-A'. 



d ^ dm 
dm " dm 



with the extended Jacobian matrix 



W (A,-i 









dm 


dAi 


dA, - 




' dm 


dm 


dm 


dAi-t 


' dm 


-El 


1 Vi, SO 



dm 



0, 



dn„, 



According to Eqs. (O and f[58j) we have = X!iLi 



— InS^'') (m) = V - — InS^'') {mi,...,mk) 
am ^ — ' ■ 



ainZ(i) 
dm 

ainz(i) 



dmi 

ainz(i) 



irii—m 
k-1 



i— 1 



9m 



ylfe-l,TO + 



OA 

ainz(i) 



Afe-i,™ 



(65) 



(66) 



(67) 



dm 



(i'.-i) 



91nZ(i) 



OA 



4— > 



o W^' Ak-2 ]o...oW' [Ai]o 



4— > 



9^ 

9m ' 



(68) 



where the last equality fohows from Eqs. ^ and ([671 . [Note that = (0, 1) 



C. Counting Optimal Packings 

In this section we attempt to determine a set of recursions to count the number of optimal packings in HNS. In 
the end, we merely succeed in providing a rigorous lower bound on the entropy density. This exercise is interesting in 
its own right as it highlights the surprising complexity in the structure of vertex covers or particle packings on this 
network. The key ingredients to provide such an approach originate with the depictions of the solutions for A: = 3 and 
4 in Figs. [7] and [51 and with the observation, in Sec.|Vl that at each finite system size iV = 2*^ + 1, exactly 3 x 2^'^ + 1 
particles can be maximally packed into the network. Let us imagine we would try to assemble the fc — 4 solutions 
from those of size fc = 3: We would have to join any two solutions at one end point and add a long link between their 
respective midpoints; the merging point becomes the new midpoint and the respective open end points remain just 
that. In the process (fc — 1) — > fc, we have to remove a single particle overall, as 

3 ^ 2('^-i)-3 _^ ij - 1 = 3 X 2'="=^ + 1. (69) 

In this construction, it appears that only the state of midpoints and end points is relevant, which we can denote by 
^nonjvnjv^ with rii G {0, 1}, depending on whether that site is (1) or is not (0) occupied by a particle. For instance, 
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Table V: Distinct classes (see the text) of solutions for HNS for each system size TV = 2* + 1. For each k, the total count adds 
up to the number of solutions given in Tab. |I] 



k 


(Oil) 


(110) 


(101) 


(111) 


3 


1 


1 


3 


2 


4 


3 


3 


10 


21 


5 


30 


30 


138 


520 


6 


4140 


4140 


22440 


162564 



the four solutions in Fig. [7]would be labeled (110) , (111) , (101) , (101) from left to right and then from top to bottom; 
we omit the reflection (Oil) of (110). In fact, a glance at Fig. [8] suggests these are the only four possibilities realized. 
We have directly enumerated these classes in Tab. |Vl 

To construct solutions of size k from those at size fc — 1, we consider all 16 pairings of these classes, which we 
symbolize by 

(^nQnN_ni±j ^nwrisN riAr^ — > (^nQnN_nN^ , (70) 

where the over-caret corresponds to the extra long-range edge added to connect the two former mid-points, prohibiting 
them from being simultaneously occupied. With that, we find the following rules: 

1. Merging two end-points into a new midpoint is possible 

(a) at no cost, when both are empty, i.e., {xxO)(Oxx)i^_^ — > (xOx)/^, making a new midpoint that is empty, or 

(b) at the expense of one particle otherwise, i.e., {xxO){lxx)f._^, {xxl){Oxx)f,_^, or {xxl){lxx)i^_-^ — > 
{xlx)^.l^ 

2. Linking the two mid-points with an edge is possible 

(a) at no cost, when at least one of the two mid-points is empty, or 

(b) at the expense of one particle, either from the left or right mid-point, if both mid-points are occupied. 

The merger can proceed only when exactly one particle gets expended, due to Eq. (|69p . Hence, the combinations of 
1(a) with 2(b) and 1(b) with 2(a) are allowed. The eight permissible mergers that are left exactly map these four 
classes onto themselves: 

[1.] (oii)rioi)fe_i ^ (oii)fe [3.] (loiKbii),-! ^ (101), [6.] (ioiKioi),_i ^ (111), 

[2.] (101)TllO)fc_i ^ (llO)fc [4.] (110Kl01),_i ^ (101), [7.] (lOlKlll),-! ^ (111), 

[5.] (110Kbll),_i ^ (101), [8.] (llfKl01),_i ^ (111), (71) 

It seems straightforward now to deduce the recursions for the number of configurations in each class, from one size 
to the next. We define the cardinality for each set as Xk = |(011),| = |(110)j,|, = |(101)j,|, and z, = |(111)^.| to 
obtain, from the rules in Eqs. STTl . 

Xk = Xk-iVk-i, (72) 
yk = 2Jk-iXk-iyk-i + '2'9k-ixl_i, 
Zk = vl-i + 2yk-iZk-i, 

with the initial conditions provided by Tab. |Vl 2:3 = l,?/3 = 3,Z3 = 2. The recursions for Xk and z, are exact, as is 
illustrated by evolving from one row to the next in Tab. IVl The recursion for y,, though, can only provide a lower 
bound on its growth. The factors of 2 in front of both terms arises from Eq. ([7T|) . as the maps [3.] and [4.] provide 
two contributions to the first while map [5.], in applying rule 2(b), gives us two ways of removing a particle in the 
second term. The "fudge factors" /, and arise because in each of these cases (and only these) the particle removal 
eliminates constraints on other particles in the respective subgraph, opening the door for an undetermined number 
of further combinations from less than optimally packed subgraphs. All we know is that these factors are larger than 
unity, but they could vary with k to an unbounded size. For further analysis, we assume that they can at least be 
approximated by constants / and g. Then, we divide the second recursion by the first in Eq. (|72p to find yk/xk ^ A 

for k — >■ 00, with X = f+ y/f'^ + 2g> 1 + v3. It is then easy to obtain asymptotically y, ~ Ax, ~ (Ais) and 
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Zk ^ 2 (Axa) (1 + 23/2/3). The total number of optimal packings is then r2fe > 2xfc + ?/fe + 2fe ^ 2^, which reduces 
to the entropy density 

..^i^>iln(A.3)>i^^^i±^« 0.1256, (73) 

using X3 = 1 and the lowest value of A. While this is a poor lower bound, it nonetheless establishes the extensivity 
of the solution-space entropy. (4^ However, its derivation also demonstrates that the structure of optimal packings is 
quite non-trivial in this network. 
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